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Abstract. Multidimensional persistence studies topological features of shapes 
by analyzing the lower level sets of vector- valued functions. The rank invariant 
completely determines the multidimensional analogue of persistent homology 
groups. We prove that multidimensional rank invariants are stable with re- 
spect to function perturbations. More precisely, we construct a distance be- 
tween rank invariants such that small changes of the function imply only small 
changes of the rank invariant. This result can be obtained by assuming the 
function to be just continuous. Multidimensional stability opens the way to a 
stable shape comparison methodology based on multidimensional persistence. 



Introduction 

The study of the topology of data is attracting more and more attention from 
the mathematical community. This challenging subject of research is motivated by 
the large amount of scientific contexts where it is required to deal with qualitative 
geometric information. Indeed, the topological approach allows us to greatly reduce 
the complexity of the data by focusing the analysis just on their relevant part. This 
research area is widely discussed in [31 H] . 

In this line of research, persistent homology has been revealed to be a key mathe- 
matical method for studying the topology of data, with applications in an increasing 
number of fields, ranging from shape description (e.g., [6l El [22l [27] ) to data simpli- 
fication [13] and hole detection in sensor networks [H] . Recent surveys on the topic 
include [T31 [Tni HHj- Persistent homology describes topological events occurring 
through the filtration of a topological space X. Filtrations are usually expressed by 
real functions ip : X —^ R. The main idea underlying this approach is that the most 
important piece of information enclosed in geometrical data is usually the one that 
is "persisienf with respect to small changes of the function defining the filtration. 

The analysis of persistent topological events in the lower level sets of the functions 
(e.g., creation, merging, cancellation of connected components, tunnels, voids) is 
important for capturing a global description of the data under study. These events 
can be encoded in the form of a parameterized version of the Betti numbers, called 
a rank invariant [5] (already known in the literature as a size function for the 0th 
homology [TSmilllT]). 

Until recently, research on persistence has mainly focused on the use of scalar 
functions for describing filtrations. The extent to which this theory can be general- 
ized to a situation in which two or more functions characterize the data is currently 
under investigation [U |5] . This generalization to vector- valued functions is usually 
known as the Multidimensional Persistence Theory, where the adjective multidi- 
mensional refers to the fact that filtrating functions are vector- valued, and has no 
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connections with the dimensionality of the space under study. The use of vector- 
valued filtrating functions, introduced in [18j . would enable the analysis of richer 
data structures. 

One of the most important open questions in current research about multidi- 
mensional persistent homology concerns the stability problem. In plain words, we 
need to determine how the computation of invariants in this theory is affected by 
the unavoidable presence of noise and approximation errors. Indeed, it is clear that 
any data acquisition is subject to perturbation and, if persistent homology were not 
stable, then distinct computational investigations of the same object could produce 
completely different results. Obviously, this would make it impossible to use such 
a mathematical setting in real applications. 

In this paper we succeed in solving this problem and giving a positive answer 
about the stability of multidimensional persistent homology (Theorem [T3|). More 
precisely, we prove that the rank invariants of nearby vector- valued filtrating func- 
tions are "close to each other" in the sense expressed by a suitable matching dis- 
tance. We point out that, until now, the stability of persistent homology has been 
studied only for scalar-valued filtrating functions. 

Our stability result for the multidimensional setting requires us to use some 
recently developed ideas to investigate Multidimensional Size Theory p] . The proof 
of multidimensional stability is obtained by reduction to the one-dimensional case, 
via an appropriate foliation in half-planes of the domain of the multidimensional 
rank invariants, and the definition of a family of suitable and possibly non-tame 
continuous scalar-valued filtrating functions. Our result follows by proving the 
stability property in this one-dimensional case. We observe that the results obtained 
in [9], for tame scalar functions, and in [51, under some finiteness assumptions, 
cannot be applied here. Indeed, we underline that our stability result requires the 
functions only to be continuous. Our generalization of one-dimensional stability 
from tame to continuous functions is a positive answer to the question risen in [5] . 

Our main point is that all the proofs carried out in |17) to analyze size functions 
via diagrams of points, the so called persistence diagrams, as well as those in p^fTl] 
to prove the stability of size functions, can be developed in a completely parallel 
way for rank invariants associated with continuous functions. As a consequence, 
we do not repeat the technical details of all the proofs, whenever the arguments 
are completely analogous to those used in the literature about size functions. For 
the same reason, we refer the reader to [I] (see also |5) for the proofs of the 
Reduction Theorem 13.21 and Theorem [33] (Stability w.r.t. Function Perturbations) 
used to deduce our main result about the stability of multidimensional persistent 
homology. 

Finally, we emphasize that throughout this paper we work with Cech homology 
over a field. In the framework of persistence, Cech homology has already been con- 
sidered by Robbins in [211 [25]. In our opinion, a strong motivation for using Cech 
homology in studying persistent homology groups is that, having the continuity 
axiom, it ensures that the rank invariant can be completely described by a persis- 
tence diagram, unlike the singular and simplicial theories that guarantee that such 
a description is complete only outside a set of vanishing measure, as explained in 
Section [2. 2. II We point out that the Cech approach to homology theory is currently 
being investigated for computational purposes [23] . 

1. Basic Definitions and the Main Result 

In this paper, each considered space is assumed to be triangulable, i.e. there 
is a finite simplicial complex with homeomorphic underlying space. In particular, 
triangulable spaces are always compact and metrizable. 
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The following relations :< and -< are defined in M": for u = (ui, . . . and 
V = . . . , Vn), we say u ^ v (resp. u ~< v) ii and only if Ui < Vi (resp. Ui < Vi) 
for every index i = 1, . . . ,rt. Moreover, R" is endowed with the usual max-norm: 
||(mi,U2, ■ • ■ ,Un)||oo = maxi<i<„ \ui\. 

We shall use the following notations: A+ will be the open set {{u, v) G M." x M" : 
u -< v}. Given a triangulable space X, for every n-tuple u = {ui, . . . , u„) G M" and 
for every function (p : X ^ R", we shall denote by X{(f ^ u) the set {x ^ X : 
ipi{x) <Ui, i = I, . . . ,n}. 

The definition below extends the concept of the persistent homology group to a 
multidimensional setting. 

Definition 1.1. Let tt^"'"'' : Hk{X{(p ^ u)) Hk{X{0 < v)) be the homomor- 
phism induced by the inclusion map tt^"'"^ : X{ip ^ u) ^ X{0 ^ v) with u ^ v, 
where Hk denotes the fcth Cech homology group, li u < v, the image of 7r["'"^ is 
called the multidimensional kth persistent homology group of {X,ip) at {u,v), and 
is denoted by H^'^\x,Lp). 

In other words, the group iJ^"'"-' (X, i^) contains all and only the homology classes 
of cycles born before u and still alive at v. 

For details about Cech homology, the reader can refer to [E] . 

In what follows, we shall work with coefficients in a field K, so that homology 
groups are vector spaces, and hence torsion-free. Therefore, they can be completely 
described by their rank, leading to the following definition (cf. [5]). 

Definition 1.2 (feth rank invariant). Let X be a triangulable space, and : X ~* 
R" a continuous function. Let fc G Z. The kth rank invariant of the pair (X, ip) 
over a field IK is the function p(x,0),k ■ — > N U {oo} defined as 

P(x,y),fc(",'y) = rankTT^"'"''. 

By the rank of a homomorphism we mean the dimension of its image. We shall 
prove in Lemma 1 2 . 1 1 that . under our assumptions on X and ip, the value oo is never 
attained. 

The main goal of this paper is to prove the following result. 

Theorem 1.3 (Multidimensional Stability Theorem). Let X be a triangulable 
space. For every fc G Z, there exists a distance Dmatch on the set {p(x,ip),k I <P '■ 
X M" continuous} such that 

(^PiX,(p).,k, P(x.tf),k^ — ^^ll'^(^) ^ '4'{^)\\oo■ 
^he construction of Dmatch will be given in the course of the proof of the theorem. 

1.1. Idea of the proof. Since the proof of Theorem 11.31 is quite technical, for 
the reader's convenience, we now summarize it in its essential ideas. We inform 
the reader that most of the intermediate results needed to prove Theorem ll.3| 
as well as the overall fiow of the proof, perfectly match well-established results 
developed in Size Theory. For these reasons, whenever this happens, details will 
be skipped, providing the reader with suitable references. However, in the present 
paper, terminology will stick to that of Persistence Theory as much as possible. 

The key idea is that a foliation in half-planes of A+ can be given, such that the 
restriction of the multidimensional rank invariant to these half-planes turns out to 
be a one-dimensional rank invariant in two scalar variables. This approach implies 
that the comparison of two multidimensional rank invariants can be performed leaf 
by leaf by measuring the distance of appropriate one-dimensional rank invariants. 
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Therefore the stabihty of multidimensional persistence is a consequence of the one- 
dimensional persistence stability. 

As is well known, in the case of tame functions, the one-dimensional persistence 
stability is obtained by considering the bottleneck distance between persistence 
diagrams, i.e. finite collections of points (with multiplicity) lying in above the 
diagonal [S]. We show that the same approach works also for continuous functions, 
not necessarily tame. We recall that a filtrating function is said to be a tame 
function if it has a finite number of homological critical values, and the homology 
groups of all the lower level sets are finitely generated. The main technical problem 
in working with continuous functions, instead of tame functions, is that persistence 
diagrams may have infinitely many points, accumulating onto the diagonal. This 
difficulty is overcome following the same arguments used in [llj for proving stability 
of size functions. 

2. Basic properties of the rank invariant and one-dimensional 
stability for continuous functions 

The main result of this section is the stability of the one-dimensional rank in- 
variant for continuous functions (Theorem 12. 19p . It generalizes the main theorem 
in [S] that requires continuous tame functions. Its proof relies on a number of basic 
simple properties of rank invariants that are completely analogous to those used in 
[TT], fT7] to prove the stability of size functions. 

2.1. Properties of the multidimensional rank invariant. The next Lemma 
12.11 ensures that the multidimensional fcth rank invariant (Definition II. 2p is finite 
even dropping the tameness condition requested in [9] . We underline that the rank 
invariant finiteness is not obvious from the triangulability of the space. Indeed, the 
lower level sets with respect to a continuous function are not necessarily triangulable 
spaces. 

Lemma 2.1 (Finiteness). Let X be a triangulable space, and (p : X —>■ M" a 
continuous function. Then, for every {u,v) G A+, it holds that P{x,ip),k{u,v) < 

+ (X). 

Proof. Since X is triangulable, we can assume that it is the support of a simplicial 
complex K and that a distance d is defined on X . 

Let us fix (m, w) £ A+, and choose a real number e > such that, setting 
e — {e,...,e),ii+2e~<v. 

We now show that there exist a function ip : X ^ K", a subdivision K" of K, 
and a triangulation L" of ip{X), such that (i) the triple {ijj, K" , L") is simplicial, 

and (a) maxj^gx - ■>pix)\\ao < £■ 

Indeed, by the uniform continuity of each component ipi of (p, there exists a real 
number 6 > such that, for i — 1, . . . , n, \(pi{x) — ipi{x')\ < e, for every x,x' £ X 
with d{x,x') < S. We take a subdivision K' of K such that mesh(i4r') < 6, and 
define tp{^) = ^i^) for every vertex x of K' . Next, we consider the linear extension 
of Tp to the other simplices of K' . In this way, V' is linear on each simplex of K' . 

Since ip is piecewise linear, ip{X) is the underlying space of a simplicial complex 
L' . By taking suitable subdivisions K" of K' and L" of L', tp also sends simplices 
into simplices and therefore {ip, K" , L") is simplicial (cf. Thm. 2.14]). This 
proves (z). 

To see (m), let us consider a point x belonging to a simplex in K', of vertices 
vi, . . . ,Vr. Since x = ' '^i' with Ai, . . . , > and J2i=i — 1' ^^'^ V' 

is linear on each simplex, it follows that [[(^(x) — V'(2^)||oo — ~ ^l=i^i ' 
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V'K)|loo = IM(^) " ELi^* • ^("i)|L = II ELi^i • 'fix) - • ^K)IL - 

ELi •^ill'^(a^) - <^(«^)||oo < 

We now prove that, since K" , L") is simphcial, it holds that Hk[X{tp < u+e)) 
is finitely generated. Indeed, since the intersection between a simplex and a half- 
space is triangulable, there exists a subdivision L'" of L" such that ip{X) Ci {x G 
K" : X ^ u + e} is triangulated by a subcomplcx of L'". By jlHl Lemma 2.16], there 
is a subdivision X'" of such that {tp, K'" , L'") is simplicial. It follows that 
X{tp <u + e) is triangulable, and hence Hk{X{ip d: u + e)) is finitely generated. 

Since u + 2e -< v and maxx^x — 4'{x)\\oo < e, we have the inclusions 

X{0 < u) ^ X{%Ij < u + ^ ^ X{0 < v) , inducing the homomorphisms Hk{X{0 < 

u)) ^ Hk{X{-if d u + e)) ^ Hk{X{Lp< v)). By recalling that Hk{X{tP ^ u + e)) 
is finitely generated, and since rank(jfe o ifS) < rank(jfc), we obtain the claim. □ 

The following Lemma 12.21 and Lemma 12.31 generalize to the multidimensional 
setting analogous results valid for n — 1. We omit the trivial proof of Lemma 12.21 

Lemma 2.2 (Monotonicity). P[x.0).k{''^i'^) is non- decreasing in the variable u and 
non-increasing in the variable v. 

Lemma 2.3 (Diagonal Jump). Let X, Y be two triangulable spaces, and f : X ^ Y 
a homeomorphism. Let (p : X —f M", : Y ^ M" be continuous Junctions such 
that maxa;gx Hv'l^;) — tp o /(a;)||oo < h. Then, setting h — {h, . . . ,h), for every 
{u,v) e A+, we have P{x,(p),k{u- h,v + h) < w). 

Proof. Since \\ip — ip o /||oo < h, we have the following commutative diagram 
Hk{X{0 ^ u - h)) Hk{X{0 <v + h)) 



Hk{Y{iP < u)) '-^-^ Hk{Y{iP d V)) 

where ik and jk are induced by inclusions, and the vertical homomorphisms are 
induced by restrictions of / and /"^, respectively. Thus the claim follows by ob- 
serving that rank ik < rank j^. □ 

2.2. Properties of the one-dimensional rank invariant. Now we confine our- 
selves to the case n — 1. Therefore, for the sake of simplicity, the symbols ip,u,v 
will be replaced by (p,u,v, respectively. We remark that A+ reduces to be the 
set {{u,v) e : u < v}. Moreover, we use the following notations: A = 9A+, 
A* = A+ U {(^,00) : u g M}, and A* = A* U A. Finally, we write \\ip\\oa for 
maxj,gx \vix)\. 

2.2.1. Right- continuity of the one- dimensional rank invariant. In what follows we 
shall prove that, using Cech homology, the one-dimensional rank invariant is right- 
continuous with respect to both its variables u and v, i.e. lim„^i(+ P(x,ip).k{u,v) = 
P(x,¥'),fc('^, and lim^-,v+ P(x,ip),k{u,v) = P(x,v),k{u,v). This property wiU be 
necessary to completely characterize the rank invariant by a persistence diagram, 
a descriptor whose definition will be recalled later in this section. In the absence of 
right-continuity, persistence diagrams describe rank invariants only almost every- 
where, thus justifying the use of Cech homology in this context. 

The next example shows that the right-continuity in the variable u does not 
always hold when persistent homology groups are defined using simplicial or singular 
homology, even under the tameness assumption. 
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Figure 1. A lower level set X{lp < u), for a sufficiently small 
u > 0, as considered in Example 12.41 corresponds to a dilation 
(shaded) of our Warsaw circle. 



Example 2.4. Let X be a closed rectangle of containing a Warsaw circle (see 
Figure [T]). Let also (p : X ^M.he the Euclidean distance from the Warsaw circle. 

It is easy to see that (p is tame on X (with respect to both singular and Cech ho- 
mology) . Moreover, the rank of the singular persistent homology group h["'"^ {X, ip) 
is equal to 1 for v > u > and v sufficiently small, whereas it is equal to when 
w = 0, showing that singular persistent homology is not right-continuous in the 
variable u. 

Analogously, it is possible to show that simplicial or singular homology do not 
ensure the right-continuity in the variable v (see Appendix [A]) . 

Let us fix two real numbers u < v and, for u < u' < u" < v, consider the 
following commutative diagram 



Hk{X{^ < u')) ^ Hk{X{>f < u")) 

Hk{X{ip < v)) ^ Hki^iv < ^))- 



By recalling that H^^'^Xx.Lp) — im TT^"'" , from the above diagram ([T]), it is easy 
to see that each tt^" induces the inclusion map a^" : (X, Lp) —>■ 

jji^u '^)(^x^ip). The following Lemma [2.51 states that, for every u" > u' > u, with 
u" sufficiently close to u, the maps ct^" are all isomorphisms. 

Lemma 2.5. Let {u,v) G A+, and let 4"'^""' : Hi""' '''\x, ip) iJ^""'"^ (X, 93) 

be the inclusion of vector spaces induced by the map . Then there exists it, 

with u < u < V, such that the maps aj!^ are isomorphisms for every u' , u" with 
u < u' < u" < u. 

Proof. We observe that the maps cr^" ■* are injective by construction (see diagram 



([!])). Moreover, by the Finiteness Lemma [2.11 and the Monotonicity Lemma | 
there exists u, with u < u < v, such that P{x,ip).kW t^) '■— rank h'^ '^■'(X, 1^9) is 
finite and equal to P(x,<p),kW , v) := rank ''"\x, if) whenever u < u' < u" < ii. 
This implies that cr^" are isomorphisms. □ 
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Analogously, by considering the commutative diagram 

Hk{X{^ < u)) Hk{X{^ < u)) 

Jii..v') 

we obtain induced maps r^^ ^ : Hj^'^ \x,(p) -ff^"'" \x,(p), and prove that 
they are isomorphisms whenever v' , v" are sufficiently close to v, with v < v' < v" . 

Lemma 2.6. Let {u,v) e A+, and let r^"''""^ : Hj^''''\x,(p) Hj^'''"\x, ip) be 
the homomorphism of vector spaces induced by the map tt^ ^ . Then there exists 

V > V such that the homomorphisms t^" ■* are isomorphisms for every v' , v" with 

V < v' < v" < V . 



Proof. The proof is essentially the same as that of Lemma 12.51 after observing 
that the maps r^" ^ are surjections between vector spaces of the same finite 
dimension. □ 

Lemma 2.7 (Right-Continuity). p(x.ip),ki'^T^) right- continuous with respect to 
both the variables u and v. 

Proof. In order to prove that lim,(_,ji+ P(x,ip).k{u,v) = P(x.ip),k{^,v), by the Mono- 

tonicity Lemma [521 it will suffice to show that Hj^^'^\x, ip) = H^'^\x, ip), where 
it is taken as in Lemma 12.51 To this end, we consider the following sequence of 
isomorphisms 

Hj. (A, 93) = im TTj. ' = im hmTTj. 

^ limim 4"'''^) = lim7?^"''"^(X,(^) i/f ^'^^X,^). 

Let us now show how these equivalences can be obtained. 

By the continuity of Cech Theory (cf. ^ Thm. X, 3.1]), it holds that im 4"'"'' = 
im limTT^" where limTr^" is the inverse limit of the inverse system of ho- 
momorphisms 7r[" : flk{X{(p < u')) Hk{X{ip < v)) over the directed set 
{u' € M. : u < u' < 11} decreasingly ordered. 

Moreover, since the inverse limit of vector spaces is an exact functor, it pre- 
serves epimorphisms and hence images. Therefore, it holds that im limTr^" '"^ = 

limim 7r[" — limij[" ''"\x, ip), where the last inverse limit is taken with respect 

to the inverse system (^Hj^!^ '^\x,p),a^ ■'^ over the directed set {u' G M : w < 

u' < u} decreasingly ordered, and ct^" are the maps introduced in Lemma 12.51 

Finally, lim7J^"''"^(X, ip) ^ h]^''^'\x, ip). Indeed, hmiJ^"''"^(X, ip) is the inverse 

limit of a system of isomorphic vector spaces by Lemma 12.51 

Analogously for the variable w, applying Lemma l2.6l □ 



2.2.2. Stability of the one- dimensional rank invariant. The following Lemma 
and Lemma 12.91 can be proved in the same way as the analogous results holding 
when A: = (see [TTj). 

Lemma 2.8. The following statements hold: 
(i) For every u < m\\\ip, P{x,<p),k{u,v) = 0. 
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(a) For every v > Taaxip, P(x,ip),k{''^Tv) is equal to the number of classes in 
Hk{X), having at least one representative in X{(p < u). 

Since, for ui < U2 < vi < V2, the number of homology classes born between ui 
and U2 and still alive at vi is certainly not smaller than those still alive at V2, we 
have the next result. 

Lemma 2.9 (Jump Monotonicity). Let ui,U2,vi,V2 be real numbers such that 
Ui < U2 < Vi < V2. It holds that 

P(X,ip)Au2, Vl) - P(X,ip),fe(Ml, ""l) > P(X,v3),fc(M2, V2) - P(X.,ip).k{ui,'V2)- 

Lemma lT^ iustifies the following definitions of multiplicity. Since we are working 
with continuous instead of tame functions, we adopt the definitions introduced in 
[T7j rather than those of [9J . Although based on the same idea, the difference relies 
in the computation of multiplicity on a varying grid, instead of a fixed one. So 
we can work with an infinite number of (possibly accumulating) points of positive 
multiplicity. Due to the lack of a well-established terminology for points with 
positive multiplicity, we call them cornerpoints, which is normally used for size 
functions. 

Definition 2.10 (Proper corncrpoint). For every point p = {u, v) E A+, we define 
the number p,k{p) as the minimmn over all the positive real numbers e, with u-\-e < 
V — e, oi 

P(x,^),k{u+e, v-e)-p(^x..^).k{u-e,v-e)-p(x,^),k{'^+£, v+e)+ p(x,^),k{u-e, v+e). 

The number p,k{p) will be called the multiplicity of p for P(^x,ip),k- Moreover, we 
shall call a proper cornerpoint for P(x.ip),k any point p e A"*" such that the number 
Mfe(p) is strictly positive. 

Definition 2.11 (Cornerpoint at infinity). For every vertical line r, with equation 
w = u, M G M, let us identify r with (m, 00) S A*, and define the number Hkij") as 
the minimum over all the positive real numbers e, with u + e < 1/e, of 

P(X.,v).k [u + - P{X.ip),k (^U ~ ^' ~^ ■ 

The number pkif) will be called the multiplicity of r for P(x.ip).k- When this finite 
number is strictly positive, we call r a cornerpoint at infinity for P(x.ip).k- 

The concept of cornerpoint allows us to introduce a representation of the rank 
invariant, based on the following definition [5]. 

Definition 2.12 (Persistence diagram). The persistence diagram Dk{X,(p) C A* 
is the multiset of all cornerpoints (both proper and at infinity) for P(x,ip),k, counted 
with their multiplicity, union the points of A, counted with infinite multiplicity. 

In order to show that persistence diagrams completely describe rank invariants, 
we give some technical results concerning cornerpoints. 

The Monotonicity Lemma \2.2\ and Lemmas 12.71 (Right-Continuity) and 12.91 
(Jump Monotonicity) imply the following result, by the same arguments as in |17| . 

Proposition 2.13 (Propagation of Discontinuities) . If p — {u,v) is a proper cor- 
nerpoint for P(x,ip),k: then the following statements hold: 

(i) If u < u < V, then v is a discontinuity point for P(x,ip),kiu., ■); 
(ii) If u < V <v, then u is a discontinuity point for P(x,f ),k{'i'^)- 
If f = (m, 00) is a cornerpoint at infinity for P(x.ip).kj then it holds that 
(Hi) If u < V, then u is a discontinuity point for P(x,ip),k{'i 
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We observe that any open arcwise connected neighborhood in A+ of a disconti- 
nuity point for P(x,ip).k contains at least one discontinuity point in the variable u or 
V. Moreover, as a consequence of the Jump Monotonicity Lemma l2.91 discontinuity 
points in the variable u propagate downwards, while discontinuity points in the 
variable v propagate rightwards. So, by applying the Finiteness Lemma 12.11 we 
obtain the proposition below. 

Proposition 2.14. Let A; G Z. For every point p — {u,v) G a real number 
e > exists such that the open set 

Weip) = {(m, v) : \u ~ Hi < e,\v - v\ < e,u ^ u,v ^ v} 

is contained in A'^ , and does not contain any discontinuity point for P[x,f),k- 

As a simple consequence of Lemma 12.81 and Proposition 12.131 (Propagation of 
Discontinuities), we have the following proposition. 

Proposition 2.15 (Localization of Cornerpoints). Ifp — {u, v) is a proper corner- 
point for P(x.ip).k, then p G {(w, v) G A"*" : minip < u < v < maxyj}. 

By applying Proposition 12.131 and Proposition 12.141 and recalling the finiteness 
of P(x,ip),k (Finiteness Lemma [2. ip . it is easy to prove the following result. 

Proposition 2.16 (Local Finiteness of Cornerpoints). For each strictly positive 
real number e, P(x.ip),k has, at most, a finite number of cornerpoints in {{u,v) G 
R"^ -.u + e <v}. 

We observe that it is easy to provide examples of persistence diagrams containing 
an infinite number of proper cornerpoints, accumulating onto the diagonal A. 

The following Theorem 12 . 1 71 shows that persistence diagrams uniquely determine 
one-dimensional rank invariants (the inverse also holds by definition of persistence 
diagram). It is a consequence of the definitions of multiplicity (Definitions 12.101 
and I2.11| . together with the previous results about cornerpoints, and the Right- 
Continuity Lemma 12.71 in the same way as done in [17] for size functions. We 
remark that a similar result was given in under the name of k-triangle Lemma. 
Our Representation Theorem differs from the fc-triangle Lemma in two respects. 
Firstly, our hypotheses on the function ip are weaker. Secondly, the A:-triangle 
Lemma focuses not on all the set A+, but only on the points having coordinates 
that are not homological critical values. 

Theorem 2.17 (Representation Theorem). For every (u, u) G A+, we have 
P(x,<p),k{u,v) = pk{{u,v)). 

(u,i.)eA* 

As a consequence of the Representation Theorem 12.171 any distance between 
persistence diagrams induces a distance between one- dimensional rank invariants. 
This justifies the following definition OE]. 

Definition 2.18 (Matching distance). Let X be a triangulable space endowed with 
continuous functions (p,^ : X R. The matching distance dmatch between P(^x,ip).k 
and P[x,iii),k is equal to the bottleneck distance between Dk{X,ip) and Dj^iY^ip)^ 
i.e. 

(2) rfmatch (p(x,¥'),fe,P(x,v),fc) = inf max \\p - l{p)\\^, 

where 7 ranges over all multi-bijections between Dk{X,(p) and Dk{X,t}}), and for 
every p = {u,v),q = {u' , v') in A*, 

Hp ~ 9II00 = mm < max {\u ~ u'\, \v — v'\} , max 
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with the convention about points at infinity that oo — y = y — (X) = oo when y ^ oo, 
oo — oo = 0, ^ = oo, |oo| = oo, min{c, 00} = c and max{c, 00} — 00. 



In plain words, || • II55 measures the pseudodistance between two points p and q 
as the minimum between the cost of moving one point onto the other and the cost 
of moving both points onto the diagonal, with respect to the max-norm and under 
the assumption that any two points of the diagonal have vanishing pseudodistance. 
We observe that in ^ we can write max instead of sup, because by Proposition 
12.151 (Localization of Cornerpoints) proper cornerpoints belong to a compact subset 
of the closure of A+. 

We are now ready to give the one-dimensional stability theorem for the rank 
invariant with continuous functions. The proof relies on a cone construction. The 
rationale behind this construction is to directly apply the arguments used in Hj for 
size functions, eliminating cornerpoints at infinity, whose presence would require us 
to modify all the proofs. 

This stability theorem is a different result from the one given in [3] , weakening the 
tameness requirement to continuity, and actually solving one of the open problems 
posed in that work by the authors. 

Theorem 2.19 (One-Dimensional Stability Theorem). Let X be a triangulable 
space, andip,ip : X two continuous functions. Then dmatch{p{x,ip).kT P{x,ip),k) ^ 

ll'/f- V'lloo- 

Proof. In what follows we can assume that X is connected. Indeed, if X has r 
connected components Ci, . . . , C^, then the claim can be proved by induction after 
observing that Dk{X, ip) — Dk{Ci, V5|Ci)j for every fc e Z. 
For fc = 0, the claim has been proved in [TTl Thm. 25]. 

Let us now assume fc > 0. Consider the cone on X , X = {X x I) / {X x {1}) (see 
Figure HI). 







p< = 1 






- 1 ^ 0/ 




<rx 



Figure 2. The cone construction used in the proof of Theorem 
12.191 The cycles in the cone are trivial. 



Since X is triangulable, so is X. We also consider the continuous function 
(p : X ^ M. taking the class of [x. t) to the value <p{x) • (1 — -I- M • i, where 
Af = 3 • (maxji^l + max 1-01) + 1. This choice of M, besides guaranteeing that 
(m, M) G A+ when u < max max \ip\, will be useful later. 

By construction, it holds that 



P(x,y),fc("'«) = I 0, 



if w > M. 



Indeed, it is well known that X is contractible (see [501 Lemma 21.13]), explaining 
why Pf^^ ^-j J, (it, v) — when v > M . The other part of the claim follows from the 
observation that, for every v < M, identifying X{(p < v) x {0} with X{(p < v), the 
lower level set X{(f < v) is a strong deformation retract oi X{(p < v). To see this, 
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it is sufficient to consider the obvious retraction r : (x, t) i— > cc and the deformation 
retraction H : X{ip < v) x I ^ X{ip < d), H{{x, t), s) = {x,t-{l- s)). This yields 
the following commutative diagram 

HkiX{<f < u)) Hk{X{f < u)) 



Hk{X{^ < v)) Hk{X{ip < v)), 

where the horizontal maps are isomorphisms, so that Hj^^'^\x,ip) = Hj^^'^\x,ip) 
when V < M. 

Clearly, a point of A+ is a proper cornerpoint for P(x,ip),k if and only if it is a 
proper cornerpoint for p^^ ^, with the ordinate strictly less than M. Moreover, 
a point (u, oo) of A* is a cornerpoint at infinity for P(x.ip),k if and only if the 
point {u, M) G is a proper cornerpoint for P(^x (p) k- remark that there are 
no cornerpoints (m, w) for P(^x 0) k '^iien maxjiysj < v < M. Analogously, we can 
construct ip : X R out of ip with the same properties. 

Following the same technical steps as in [TTl Thm. 25], simply substituting 
P(x,^),k for P(x,v),o> it is possible to prove the inequality d,natch{pt^x,(p),k^ P(x,^),k) ^ 
llf? — V'lloo- To this end, we need to apply Lemma [T^ and Propositions 12. 131 12.141 
12.151 and l2. 161 Therefore, since — -01100 = ||</5 ~ V'lloo, it is sufficient to show that 

dmatch{P{X,ip),k-, P(X,iP),k) < '^match{p (^x ,(p) ,k^ P (X ,k^ ■ 

We can prove that a multi-bijection 7 between Dk{X,(p) and Dk{X,ip) exists, 
with rfmatc/»(P(x,y),fc. P(x j ;),/c) = ^^^peDk{x,v) " 7(p)lls5- This can be done by 
applying Proposition (^HH as in [TTl Thm. 28]. Such a 7 will be called optimal. 

Since 7 is optimal, then 7 takes each point {u,v) e Dk{X^ip), with v = M, 
to a point {u',v') £ with v' — M. Indeed, if it were not true, i.e. 

7((m, M)) = {u',v') with v' < M, then v' < max|-0|, and by the choice of M, 
we would have ||(m,M) — 7((u,M))||55 > ||(^ — V'lloo, contradicting the fact that 
dmatch{P(x ,0),k^ P(x ,i,),k) ^ llv' - V'lloo- Analogously for 7"^ proving that 7 maps 
cornerpoints whose ordinate is smaller than M into cornerpoints whose ordinate is 
still below M. 

Let us now consider an optimal multi-bijection 7 between Dk{X, (p) and Dk{X, 
We now show that there exists a multi-bijection 7 between Dk{X, ip) and Dk{X, 
such that maxpg£,fc(x,¥>) lb - 7(p)ll55 = max^ez5fc(x,0) \\P ~ 7(p)lls5, thus proving 
that d^atch{p(x,^),k,P(x,^),k) < dmatch{P(^x,(p),k^P{x,ii),k)- Indeed, we can de- 
fine 7 : Dk{X,(p) Dk{X,^) by setting j{{u,v)) = j{{u,v)) if w < 00, and 
"f{{u,v)) = {u',v), where u' is the abscissa of the point 7((u, A/)), if w = 00. This 
concludes the proof. □ 

3. Proof of the Multidimensional Stability Theorem 11.31 

We now provide the proof of the Multidimensional Stability Theorem 11.31 It 
will be deduced following the same arguments given in [l] to prove the stability of 
multidimensional size functions. Proofs that are still valid for fc > without any 
change will be omitted. 

We start by recalling that the following parameterized family of half-planes in 
M" X R" is a fofiation of A+. 

Definition 3.1 (Admissible pairs). For every unit vector I — (Zi,...,^„) of K" 
such that > for z = 1, . . . , n, and for every vector b = (5i, . . . , 6„) of R" such 
that X]r=i ~ 0' shall say that the pair {I, b) is admissible. We shall denote the 
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set of all admissible pairs in M" x M" by Adirin. Given an admissible pair (/, 6), we 
define the half-plane Ti'y-g') of E" x M" by the following parametric equations: 

{u = si + b 
V = tl + b 

for s,t eR, with s <t. 

The key property of this foliation is that the restriction of P{x.(p),k to each leaf 
can be seen as a particular one-dimensional rank invariant, as the following theorem 
states. 

Theorem 3.2 (Reduction Theorem). Let {l,b) be an admissible pair, and F'^-^ : 
X —I- R be defined by setting 

ipiix) - i 



il,b)^ i=l,...,n [ 1^ 

Then, for every (u, v) = (si + b,tl + b) G '^(fb) following equality holds: 
P(x,y),fc(", v) = . ^(s, t) . 

^ ' (l,b)" 

As a consequence of the Reduction Theorem 13.21 we observe that the identity 

P(X,0),k = P(x th) k if O'^ly if d,natch{P(x F"^ ) k^ P I v W ) = ^^^'"^ 

admissible pair 

The next theorem gives a stability result on each leaf of the foliation. It repre- 
sents an intermediate step toward the proof of the stability of the multidimensional 
rank invariant across the whole foliation. 

Theorem 3.3 (Stability w.r.t. Function Perturbations). If X is triangulable and 
Lp,ip : X ^ R" are continuous functions, then for each admissible pair {l,b), it 
holds that 

0'm.atch\P(xpt^'\.k^P<YW^ \ l.) - 



We are now ready to deduce the Multidimensional Stability Theorem ll.3l 

Proof, (of Theorem\r^ We set D„iatch{P{x,<p),k, P(^x.^),k) = ^^P{Lb)i£Adm„ ^^^^ ' 

dmatchiPix F-^ k' P < V t?^ \ i) ' Then, Dmatch is a distance on {p(x.^),k I <^ : -'^^ ^ 

' (i.b)'' '•-^'-^(r.S)-'' 

M" continuous} that clearly proves the claim. □ 

Roughly speaking, we have thus proved that small changes in the vector-valued 
filtrating function induce small changes in the associated multidimensional rank 
invariant, with respect to the distance Dmatch- 
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A. Appendix 

The next example shows that the rank invariant is not right-continuous in the 
variable v when singular or simplicial homology are considered instead of Cech 
homology. We recall that an example concerning the right-continuity in the variable 
u has been given in Example 12.41 
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Example A.l. Let 5* C be a sphere parameterized by polar coordinates {a,f3), 
— f < a < f and P € [0, 2tt). For every (3 e [0, 27r), consider on S the paths 7^ : 
(-§ ,0) ^ S and 7| : (0, f ) 5 defined by setting, for i ^ 1,2, 7^(a) = (a',/3') 
with a' = a and /?' = (/? + cot a) mod 27r. We observe that each point of the set 
S* = {{a, P) & S : a ^ A \a\ ^ ^} belongs to the image of one and only one path 
7^. Such curves approach more and more a pole of the sphere on one side and the 
equator, winding an infinite number of times, on the other side (see, for instance, 
in Figure [3] (a) , the paths 7^ and ^'i^ lying in the northern hemisphere) . 

A+ 



u 

^(P) = ^(Q) 

(«) (6) 

Figure 3. (a) Two of the paths covering the northern hemisphere 
considered in Example lA.il (&) The 0th rank invariant of the pair 
(5*, . On the discontinuity points highlighted in bold red, the 
0th rank invariant computed using singular homology takes value 
equal to 2, while using Cech homology, the value is equal to 1, 
showing the right-continuity in the variable v. 




Then define the C°° function (^* : S** — > R that takes each point P = 7^ (a) G S* 

to the value exp ( -. — ^ — ^ ) sin(/3). Now extend (p* to a C°° function (^9 : 5 ^ M 

V "Hf-I"lj J 

in the unique possible way. In plain words, this function draws a ridge for [3 G (0, tt), 
and a valley for (3 £ (tt, 27r). Moreover, observe that the points P = ^) and 
Q = (— f , ^) of the sphere are the unique local minimum points of ip. 

Let us now consider the 0th rank invariant of the pair (5, Lp) . Its graph is 
depicted in Figure [3] (6). The points P and Q belong to the same arcwise connected 
component of the lower level set S{ip < e) for every e > 0, whereas they do not for 
e — 0, since the paths 7^ {i = 1, 2) are an "obstruction" to construct a continuous 

path from P to Q. Hence, the singular rank invariant P(^s.'^).o is not right-continuous 
in the second variable at u = 0, for any u with minp) < u < 0. 
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